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The Strong Coupling Fixed-Point Revisited 
A.C. Hewson* 
Dept. of Mathematics, Imperial College, London SW7 2BZ, UK. 

In recent work we have shown that the Fermi Uquid aspects of the strong coupling fixed 
point of the s-d and Anderson models can brought out more clearly by interpreting the 
fixed point as a renormalized Anderson model, characterized by a renormalized level id, 
resonance width, A, and interaction U, and a simple prescription for their calculation was 
given using the numerical renormalization group (NRG). These three parameters completely 
specify a renormalized perturbation theory (RPT) which leads to exact expressions for the 
low temperature behaviour. Using a combination of the two techniques, NRG to determine 
id, A, and U, and then substituting these in the RPT expressions gives a very efficient 
and accurate way of calculating the low temperature behaviour of the impurity as it avoids 
the necessity of subtracting out the conduction electron component. Here we extend this 
approach to an Anderson model in a magnetic field, so that id, A, and U become dependent 
on the magnetic field. The de-renormalization of the renormalized quasiparticles can then be 
followed as the magnetic field strength is increased. Using these running coupling constants in 
a RPT calculation we derive an expression for the low temperature conductivity for arbitrary 
magnetic field strength. 

KEYWORDS: Strong coupling fixed point, Anderson model, numerical renormalization group, 
renormalized perturbation theory, Fermi liquid. 



1. Introduction 

After Wilson's seminal numerical renormalization group (NRG) solution of the spin 1/2 
s-d model^ it was very soon recognized by Nozieres^ that the low energy strong coupling fixed 
point corresponds to Fermi liquid behaviour. Using a phenomenological description of the 
phase shift Nozieres then gave an analytic derivation of Wilson's result for the "x/7" ratio i?,^ 
R = 2, and derived the leading form of the low temperature conductivity in terms of the Kondo 
temperature.'^ The important characteristic of a Fermi liquid is the 1-1 correspondence of the 
single particle excitations with those of the non-interacting system. This correspondence does 
not apply in the case of the s-d model because the model has a constraint of a fixed occupancy 
nd of the impurity site, = 1 for S = 1/2, which implies a strong local interaction to enforce 
it. The Anderson model, however, which is equivalent to the s-d model {S = 1/2) in the local 
moment regime, is truly non-interacting when the interaction term U is set to zero, so the 
1-1 correspondence of the single particle excitations with those of the non-interacting system 
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should hold, even in the local moment limit. This implies that low energy fixed point of the 
s-d model and Anderson model should be described more naturally in terms of a renormalized 
Anderson model. The Anderson model"^ has the form, 

Ham = J2 ^ddUa + U nd,^nd,l + ^(Vfc4cfe,a + ^k^la^^a) + ^k,aclaCk^a: (1) 

f k,cr k,(T 

where is the energy of the impurity level, U the interaction at the impurity site, and Vfc the 
hybridization matrix element to a band of conduction electrons with energy e^.Whcn [7 = 
the local level broadens into a resonance, corresponding to a localized quasi-bound state, 
whose width depends on the quantity A(lo) = vr^^, — e^). It is usual to consider the 

case of a wide conduction band with a fiat density of states where A(u;) becomes independent 
of LO and can be taken constant A. 

In the wide conduction band limit the one-electron Green's function, Gdico) takes the form 

G<i(w) = K (2) 

where S(a;) is the self-energy and A is the width of the resonance at when [/ = 0. Near 
the Fermi level for small = S(0) -|-a;S'(0) -|- 0(a;^), and if this is substituted into (2) 

then for small io the denominator is of the same form as that for the non-interacting system 
with a renormalized level and resonance width A given by 

ed = ^(ed + S(0)), A = zA, (3) 

where z, the wavefunction renormalization factor, is given by ^ = 1/(1 — S'(0)), and the prime 
indicates a derivative with respect to The 1-1 correspondence is evident when one calculates 
the quasiparticle occupation number at T = 0, 

1 1 _i 

= 2 - - tan ( ^ I = ^du, (4) 

which is equal to the impurity occupation number ride for spin cj at T = from the Friedel 
sum rule.^ The corresponding quasiparticle density of states is given by 



[UJ - EdY + A^ 

It follows from Fermi liquid theory that the impurity specific heat coefficient 7, as calculated 
from these non-interacting quasiparticles, is exact as T — >^ 0, and is given by 

27r2 

7imp = -^-^^(O). (6) 

In the original NRG calculations for the Anderson model the low energy fixed point was 
analysed as a strong coupling V ^ 00 fixed point. ^ In this limit the impurity is decoupled, 
and so the analysis does not bring out the 1-1 correspondence of the single particle excitations 
with those of the original Anderson model. Recently we (Hewson, Oguri and Meyer^) have 
reanalysed the strong coupling fixed point as a J/ = fixed point with a finite V , which 



2/17 



J. Phys. Soc. Jp Full Paper 




Fig. 1. Plots of the parameters, A(A^)/A (crosses), €4{N)/ti/S. (triangles), lJpp{N)) / ttIS. (diamonds), 
Uhh{N)/TTA (stars) and Uph{N)/TTA (circles), with N, for the Anderson model with bare param- 
eters, U — 0.04, ttA — 0.03, and = —0.05. In this and subsequent figures the full lines are 
interpolations using the calculated data points. 



leads directly to the renormalized parameters and A. The NRG calculations are based on 
linear chain form for the Anderson model, with a discretized conduction electron spectrum 
and a discretization parameter A > 1. Starting with the impurity at the end of the chain, the 
Hamiltonians corresponding to a finite lengths of chain are diagonalized iteratively, adding a 
new site to the chain with each iteration. If the lowest energy single-particle Ep{N) or single- 
hole excitations Eh{N) of the interacting Anderson model for a chain with N + 2 sites can be 
described by a non- interacting renormalized Anderson model then E = Ep(N) or E = Efi{N) 
should satisfy the equation, 

^A-(A^-i)/2 _ g^(^) ^ Ai^-mv^Nfg^^^E), (7) 

where the effective parameters, ed(A^) and A(A^) = 7rV{N)'^ /D, should be independent of N. 
The function goo{uj) is the local Green's function for the first conduction electron site of the 
chain when decoupled from the impurity (V = 0) (see^ for details). 

In Figs. 1 and 2 we present results for ed{N) and A(A^), determined by these two equations 
as a function of N. In the first case shown in figure 1 we take a relatively weak coupling example 
with bare parameters, vrA = 0.03, = —0.05 and U = 0.04 (bandwidth 2D = 2 in all cases), 
which is such that U pd.mf (0) = 0.331 < 1, where pd,mf{0) is the mean field density of states at 
the Fermi level, Pd,mi{0) = ^/^(^Imf + ^^) 'with e^.mf = e(i + f^'^d.mf/2, and so does not satisfy 
the mean field (Hartree-Fock) criterion for a local moment. The renormalized parameters 
are independent of N for > 26, which confirms that these single-particle excitations on 
the lowest energy scale can indeed be described by a renormalized Anderson model. The 
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Fig. 2. Plots of the parameters, A(A^)/A (erosses), ed(7V)/7rA (triangles), iJpp{N)) / -k/S. (diamonds), 
Uhh{N)/iT/S. (stars) and lJph{N) /-k/S. (circles), with the N, for the Anderson model with bare 
parameters, U/nA = 5.0, ttA = 0.03, and ea = —0.05 . 



asymptotic value of A/A = 0.947 for large A'', and so A differs little from its bare value. The 
value of €(i = —0.0143, and is very approximately of the same order as that predicted from 
mean field theory c^/^j^if = -0.0166. 

In Fig. 2 we give the corresponding results for U = 0.15, keeping the other parameters the 
same. For this value of U mean field theory would predict the breaking of local spin symmetry 
as Upd^mii^) = 2.02 > 1, and hence this case corresponds to a system with a local moment and 
can be described by an effective s-d model. We see a considerable renormalization such that 
A/A = 0.0158, and becomes very small = 1.06 x 10~^ so the effective level is very close 
to the Fermi level, as to be expected in the almost localized limit ~ 1. Having determined 
ed and A, the impurity occupation and the specific heat coefficient 7 can be determined from 
eqs. (4) and (6). 

The renormalized quasiparticles must interact with one another, and this interaction must 
come into play as soon as two or more single particle excitations are created from the interact- 
ing ground state. If the lowest two-particle excitation from the ground state for the interacting 
system for a given N has an energy Epp{N), then we can calculate U by equating the energy 
difference Epp{N) — 2Ep{N) to that calculated by adding an local interaction term to the 
effective Anderson model for the non-interacting quasiparticles.®''' For finite N we can use 
this equation to define an A'^-dependent renormalized interaction Upp{N), 

Epp{N) - 2Ep{N) = ^(Ar)A(^-i)/2|^;_,(-l)p|V,* ^(-l)p, (8) 
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Fig. 3. The same plot as in Fig. 2 for A'' > 25 but shown over a smaller energy range. The parameters 
correspond to the almost localized (Kondo) limit; the renormalized level Cd/TrA (triangles) is very 
close to the Fermi level, and the limiting value of A{N)/A for large N (crosses) coincides with 
the limiting values of Upp{N))/iTA (diamonds), Uhh{N)/-KA and Uph{N)/'jTA (circles). 



where |'0p,i(— l)p is given by 

where ^'ooC'*^) is the derivative of 500 (w)- 

Alternatively we could consider the same procedure for a two hole excitation Ehh{N) and 
in a similar way define an iV-dependent renormalized interaction Uhh{N), or a particle- hole 
excitation Eph{N) to define a renormalized interaction Uph{N). In this latter case, as a positive 
U leads to particle- hole attraction, we use Ep{N) + Eh{N) — Eph{N) on the left-hand side of 
eq. (8). 

If these two particle excitations can be described by an effective Anderson model then 
Upp{N), Uhh{N) and Uph{N) should be independent of N and also independent of the particle 
and hole labels. The values of Upp{N), Uhh{N) and Uph{N), for the case U = 0.04 as a function 
of N are shown in Fig. 1, and those for U = 0.15 in Fig. 2, where ttA = 0.03 and = —0.04 
in both cases. The three interaction parameters can be seen to converge to a unique value for 
large N both in the weak and strong coupling cases. 

In Fig. 3 the renormalized interaction parameters in the strong coupling case for N > 25 
of the plot given in Fig. 2 are shown over a smaller energy range. It can be seen that Upp{N), 
Uhh{N) and Uph{N) not only have a unique limit but that this limit coincides with the 
asymptotic value of TrA{N) for large N, so that ttA = U; consequently there is only one 
effective parameter scale in the localized (Kondo) regime. 
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An alternative renormalization approach to the Wilson technique is the renormalized per- 
turbation theory (RPT)^ as originally developed to deal with the infinities arising in quantum 
electrodynamics (QED). While the Wilson approach is based on the progressive elimination 
the higher order excitations to obtain an effective low energy model, the renormalized per- 
turbation theory is essentially a reorganization of perturbation theory so that a perturbation 
expansion can be carried for the same model, but in terms of parameters appropriately renor- 
malized for very low energy scales. These low energy scales are where almost all observations 
are made in QED so that the renormalized parameters can be taken from experiment. As the 
RPT is not simply about the cancellation of infinities, but working with parameters appropri- 
ate to the energy scale under investigation, it can be applied quite generally. This approach 
has been developed for the Anderson model, and leads naturally to the quasiparticle descrip- 
tion.^' Here we briefly review some of the main results which will be required in the later 
sections of this paper. The renormalized parameters of the non-interacting quasiparticles are 
a renormalized level and resonance width A, as given in eq. (3) in terms of the self-energy 
and its derivative at T = a; = 0. The renormalized local interaction U is identified with the 
fully dressed irreducible 4- vertex, r-^^{uji,uj2,i^3,i^4:) at ui = i02 = uj^ = L04 = 0, 

U = z%^{0,0,0,0), (10) 

which is rescaled by z"^ so the propagator for the quasiparticles is normalized. 

In applying the renormalized perturbation expansion the Lagrangian for the Anderson 
model C\M{^dj ^) U) is rewritten in the form, 

£AM(ed, A, U) = CAuied, A, U) + £c(Ai, A2, A3). (11) 

The first term on the right hand side is simply the Lagrangian for the Anderson model 
expressed in terms of the renormalized parameters, and the second part is the remainder or 
counter term. The three parameters Ai, A2 and A3 are fully determined order by order in the 
expansion in powers of U by the requirement that they prevent overcounting and cancel off 
any terms that further renormalize the particles as these quantities have be taken to be fully 
renormalized already. Hence the three renormalized parameters id, A and U are sufHcient to 
specify the RPT precisely. 

The first order perturbation theory in U gives results for the impurity spin and charge 
susceptibilities at T = 0, 

xs = ^^Pdm^+upM), xc = mm-u~Pdm, (12) 

where pdi'jj) is given by eq. (5). These results can shown to be exact by the use of a Ward 
identity, and are equivalent to the results as first derived by Yamada^^ from an analysis of a 
perturbation expansion in powers of U. The 'x/7' ratio R is then given by i? = 1 -|- Upd{G). 
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The RPT result for the rcnormahzed self-energy to second order in U for the symmetric 
model gives the exact low temperature result for the conductivity cr(T) to second order in the 
temperature T, which is given by 



(7{T) = ao < 



3 VA 



ttA 



+ 0(r4 



(13) 



When the renormalized parameters are expressed in terms of the self-energy and the vertex 
function these results coincide with the exact expressions derived by Yamada and Yosida^^ 
from an analysis of perturbation theory to all orders U, and in the localized regime {U /ttA 
1) with the Fermi-liquid results of Nozieres.^ The corresponding result for the differential 
conductance through a quantum dot to second order in the applied voltage has been given by 
Oguri.^^ 

Higher order terms in U can be used to estimate the leading corrections to the Fermi 
liquid results. Estimates of the induced magnetization M{h) as a function of magnetic field 
H, where h = gfiBH/2 to order have been made using RPT to third order in U for the 
particle-hole symmetric model. The result to this order is 



M{h) 




1 + 4 




+ A 




+ B 




where the coefficients A and B are given by 

, 38 37r2 
A = — -—, and B 



644 
~9~ 



+ .. 

(14) 

(15) 



Though this estimate of the coefficient of the term is not exact, it is asymptotically 
exact for small U {U 0) and differs at the most by 4% from the Bethe ansatz result in the 
localized Kondo regime {U » ttA). 

In principle it should be possible to derive results appropriate for all energy scales for 
the RPT given the three renormalized parameters, as nothing is neglected when the counter 
terms are taken into account. However, higher order calculations get progressively more diffi- 
cult and it is unlikely that the summation of any subset of terms is likely to provide reasonably 
accurate results on all energy scales. As we shall discuss later, however, it might prove pos- 
sible to calculate a set of running coupling constants, appropriate to the energy scale under 
consideration, and use the low order RPT to cover all energy scales in this way. 

3. Renormalized Anderson Model 



The first term on the right hand side of eq. (11) can be identified as the fixed point of 
the Wilson NRG approach, because as T ^ the effect of the counter terms is to normal 
order the interaction term so that it only comes into play when two or more excitations are 
created from the ground state. The Hamiltonian for the renormalized Anderson model, which 
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Fig. 4. A plot of the renormalized parameters id, A and t/ as a function of the bare interaction U 
for an Anderson model with ttA = 0.03 and = —0.04. Over this range of U the system moves 
from a full orbital regime €d + U ^0, through an intermediate valence regime + ?7 ~ 0, to a 
localized (Kondo) regime + ?7 3> 0. Also shown are the impurity occupation number rid and the 
'x/7' or Wilson ratio R. 



describes the behaviour near the low energy fixed point, therefore, can be written as 

^AM = ^ eddUa + U : nd,^nd,i : + ^(Vfc4cfc,<7 + ^k4,ada) + efc,a4,aCfc,a> (16) 

where the colon brackets indicate that the expression within them must be normal-ordered. 
This renormalized model is similar to that used in earlier phenomenological local Fermi-liquid 
theories, but here it also includes a quasiparticle interaction term. 

We can identify U as defined from eq. (10) with that from the NRG calculation in eq. (8). 
The NRG results for id, A and U can then be substituted in eqs. (12) to evaluate the T = 
spin and charge susceptibilities. This is considerably simpler than the method used originally 
to evaluate these and is very accurate as it does not involve subtracting out the conduction 
electron component. It also follows analytically from these results that in the Kondo regime, 
Xc(0) — > and — > 1 that — > and U = vrA, as noted earlier in the results shown in Fig. 
3. It is equivalent to Nozieres' argument^ and gives the Wilson ratio R = 2. If we define the 
Kondo temperature Tk via Tk = (gfiB)'^ /^Xs{0) then U = vrA = 4Tk- 

The advantage of the renormalized Anderson approach to describe the low energy be- 
haviour of non-degenerate magnetic impurities is that all the parameter regimes; weak cou- 
pling, mixed valence, strong correlation, empty and full orbital regimes, can be described 
precisely within a single framework with, at the most, three renormalized parameters. In Fig. 
4 we take values of U over the range < U < SvrA, with the same values of and A as 
used earlier. Over this range we move from the full orbital regime for small U, through an 
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intermediate valence regime for + ?7 ~ and the Kondo regime for U ^ |ed[- We see that 
id increases at first approximately linearly with U until U/ttA 2.0, remaining very close 
to the Fermi level in the Kondo regime at higher values of U. The renormalized resonance 
width A decreases over the same range monotonically approaching zero in the limit [/ — > oo. 
The quasiparticle interaction U increases at first linearly with U, reaching a maximum for 
U/itA ~ 1.5, and then decreases so that its energy scale merges with that for A in the Kondo 
regime. The initial increase of with U in the full orbital regime is understandable in terms 
of the mean field or Hartrcc-Fock theory, which is approximately valid in this regime. In the 
mean field theory there is no oj-dependence in the self-energy so z = 1 and Amf = A. There 
is an effective level given by e^^mf = + f^'^d,mf/2, where is the total impurity occu- 

pation calculated within the mean field theory. The value of U in proportional to U for very 
small U, and for larger values of U, its value can be estimated from (8) using perturbation 
theory. These expressions give the general trend as a function of U in the full orbital regime 
Upd,mf ^ 1 but are not valid as the strong correlation regime is approached, where A begins 
to differ from the bare value A, and breaks down completely in the Kondo regime where A is 
strongly renormalized (sec reference'' for more extensive results). 

Also plotted in Fig. 4 is the Wilson ratio, which is given by i? = 1 + Upfi{0), and the total 
impurity occupation number n^. The occupation number decreases from 1.88 for [/ = to 
almost 1 in the large U regime, corresponding to localization of the d-electron. The Wilson 
ratio increases from 1 in the small U regime and asymptotically approaches 2 in the Kondo 
regime. 

4. De-renormalization as a Function of Magnetic Field 



If wc introduce a magnetic field H, wc can again generalize the definition of the renormal- 
ized parameters, e^a A and U, such that they become functions of the magnetic field H. For 
simplicity wc confine the discussion to the particlc-hole symmetric model with = —U/2. 
If we absorb the zero field Hartree-Fock contribution to the self-energy, S(0, 0) = U/2 then 
€d = 0, and S|(0, h) = — Sj^(0, h) = — S|(0, —h), where h = gp-QH/2. Eq. (3) gets replaced by 

~€A,a{h) = z{h){-ha + E^(0, h)), A{h) = z{h)A, (17) 

where z{h) is given by z{h) = 1/(1 — S'|(0, /i)). The occupation of the impurity level is still 
given by the Priedel sum rule, 

The definition (8) of U can be straightforwardly generalized to define an field dependent local 
interaction U{h). As ed^{h) = —edi{h), it will be convenient to define a single effective level 
ed{h) via ed{h) = —aidaih). The impurity magnetization M{h) at T = is then given simply 
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ln(h/7rA) 

Fig. 5. The renormalized parameters ed{h)/h. A{h)/A, [/(/i)/7rA and the Wilson ratio R{h) plotted 
as a function of the logarithm of the magnetic field h/nA for the symmetric Anderson model 
with U/ttA — 5, ttA — 0.03, and a bandwidth 2D = 2. The parameters vary from the strong 
correlation values in weak field, 7rA(0) — f/(0), R{0) = 2 = lim/,^o ^d{h)/h, to those of the bare 
model, U{h) ^ U, A{h) = A, R{h) = 1 = ed{h)/h in fields /i > [/. 



by 

M,„ _ i^,..-^ (II) ^ 

The expressions given earlier for the quasiparticle density of states (5), impurity spin and 
charge susceptibihties (12), and specific heat coefficient (6) can be generahzed to include the 
magnetic field dependence by replacing e^, A and U, by edih), A(/i) and U{h). 

The simplest way to calculate the field dependent renormalized parameters using the NRG 
for the particle-hole symmetric model is to exploit the spin-isospin symmetry of the model, 
which is such that the symmetric model in a magnetic field with positive-C/ is equivalent to 
a negative-C/ asymmetric model in the absence of a magnetic field with + U/2 = —h and 
spin and charge (isospin) interchanged. Calculations can then be carried out in the absence of 
a magnetic field using the asymmetric model with negative-C/. The isospin symmetry of the 
particle-hole symmetric model is then effectively exploited via the spin symmetry. This has 
two advantages, it enables one to retain more states in the NRG iterations, and requires no 
modification of the standard NRG program. 

In Fig. 5 we give the results for erf(/i)//i, A(/i)/A and U{h)/TrA as a function of h for 
the symmetric model with U/nA = 5, which corresponds to a Kondo temperature Tk/vtA = 
0.00204. We can follow the gradual de-renormalization of the quasiparticles as the magnetic 
field increases. In weak field id/h — > 2 as /i — > 0, so the level splitting in a magnetic field of 
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the local spin-up and spin-down states is twice the Zeeman splitting for free quasiparticles. 
This enhancement factor of 2 is the same as that for the Wilson or xll ratio. As the magnetic 
field strength increases e^{K)lh increases to a maximum, and only for magnetic field strengths 
larger than h ^ U does this ratio finally approach the free particle value, ed{h)/h — > 1. 
There is an even more dramatic increase in the value of U{h)/TTA. For h ^ we are in the 
strong correlation or Kondo limit, U{0) /it A{0) = 4TK/'7rA, so the value of U{h) is small in 
the weak field regime. However, it increases significantly with an increase of field strength to 
a maximum, approximately three times greater than the bare value of 5, before eventually 
approaching the bare value for h ^ U. 

The fact that the quasiparticle interaction becomes very large does not imply that the 
effects of the interaction become stronger; the contrary is the case. A more significant measure 
of the effect of the interactions is the product of U (h) with the quasiparticle density of states 
at the Fermi level pdiO,h). The increase in U{h) is more than compensated by the fall off 
of pd{0, h) with h as the level moves away from the Fermi level. This can be seen from the 
Wilson ratio R{h) = 1 -|- pd{h)U{h), which is also shown in Fig. 5. The value of R decreases 
monotonically from the strong correlation value 2 to the free particle value 1 in the extreme 
large field limit h U, which implies that the product pd{h)U{h) is always less than 1. The 
enhanced value of U{h) in the strong magnetic field limit, and the decrease as h ^ oo, is 
understandable in terms of mean field theory, where from eq. (8) U (h) is given by 

1 - Upd,m{{0, h) 
as P(i,mf (0, h) decreases with increase of h in this regime. 

In Fig. 6 we plot the impurity magnetization as function of the logarithm of h using 
the non-interacting quasiparticle expression in eq. (19) using the parameters ed{h) and A(/i) 
given in Fig. 5, together with the corresponding values for an s-d model with the same Kondo 
temperature. -"^^ Also plotted in the same Fig. is R{h)/4. One sees that there is complete 
agreement with the results of the s-d model until charge fluctuations begin to play a role for 
/i ~ ttA (ln(/i/7rA) ~ -3). Up to this point R{h) ^ 2 {U{h)pd{0,h) ^ 1) independent of h, 
as is well known for the s-d model,^^' but then crosses over eventually to the free electron 
value R{h) = 1 for /i » ttA. The magnetization can be seen to approach saturation more 
rapidly with h than predicted from the s-d model, due to the additional eflFect of the charge 
fluctuations. 

The calculation of the magnetization from eq. (19) does not depend on the quasiparticle 
interaction U{h), but the formula for the susceptibility in eq. (12) as a function of h does. 
A check on the values of U{h) can be made by calculating the susceptibility from eq. (12) 
using ed{h), A{h) and U{h), and comparing the result with that calculated by numerically 
diflFerentiating the results for M{h). When this check is made complete consistency is found, 
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ln(h/;tA) 



Fig. 6. The induced impurity magnetization M{h) (stars and dashed curve) as a function of ln(/i/7rA) 
for the symmetric Anderson model with the same parameter set as used in Fig. 5 {U/nA = 5) 
compared with the same values from the Bcthe ansatz solution^^' of the s-d model (full curve). 
Also shown is the Wilson ratio R{h) divided by 4 (crosses and dotted line). 



confirming the calculated values of U {h) . 

5. Perturbation Expansion with Running Coupling Constants 

The values of ed{h), A(/i) and U{h) obtained in the previous section provide a set of 
running coupling constants for a renormalized perturbation expansion for the symmetric An- 
derson model in a magnetic field. Instead of working with A and U of the model in the 
absence of a magnetic field, the perturbation expansion can be carried out using parameters 
appropriate to the field strength h. If the calculations could be carried out exactly then it 
would not matter which set of coupling constants is used, as the model is completely defined 
by any one set. However, for approximate calculations in the low energy regime the best set 
in the presence of a magnetic field h must be the set ed{h), A(/i) and U{h), because low order 
calculations in terms of these parameters give asymptotically exact results as T — >^ 0. 

We apply this approach to the calculation of the low temperature conductivity in the 
presence of a magnetic field for the particle-hole symmetric model a{T, h). The calculation is 
along the same lines as that used to derive the /i = result in eq. (13).^ All the terms in the 
RPT expansion to order U{h)'^ are taken into account, including the counter terms to this 
order. The result, when summed over the two spin components, can be written in the form, 

-f'^' = iJ^) {' + + 0(m^Krt} , (21) 
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where m{h) = 'KM{h)/ and the coefficient a2{h) is given by 




= J—^ — {1 + 2[R{h) - lYcos^m{h)} . 



(22) 



This result is exact to order if either h or U is set to zero, and we conjecture, that it is exact 
for the general case U and /t 7^ 0. In Fig. 7 we plot the logarithm of the coefficient (72 (/i) 
of the term as a function of h over a range from strongly correlated quasiparticles in weak 
field to bare quasiparticles in fields h':$> U, using the renormalized parameters shown in Fig. 
6 for the model with U/ttA = 5. In Fig. 8 we plot (72 (/i) for a more realistically obtainable 
range of h, < h < I.STk, for the same model where Tk/ttA = 0.00204. There is a very 
significant decrease in this coefficient over this magnetic field range. 

This RPT calculation has had to rely on the renormalized parameters derived from the 
NRG energy levels. It might, however, be possible to calculate them directly from eqs. (3) 
and (10), using an iterative RPT. In the asymptotically large field limit, the parameters are 
those of the bare model, but the interaction effects are small because the impurity is almost 
completely polarized. It should be possible, therefore, to use perturbation theory for the bare 
model to calculate the small renormalization effects for large h, and the set of renormalized 
parameters for this regime. With these new set of renormalized parameters, the RPT could be 
then used to calculate the renormalized parameters for slightly smaller fields, and enabling one 
to set up flow equations to continue the process to the small field regime. The perturbation 
effects at each stage should be small, as the effects of only small changes in the magnetic 
field will be calculated at each stage, and there are only three parameters, id, A and U, to 
consider, as these three fully specify the expansion. 

It might also be possible to do something similar as a function of temperature starting 
with the high temperature limit, progressively lowering the temperature, hence having run- 
ning parameters, €d{T), A(T) and U{T). We have explored this possibility by taking the 
parameters, ed{N), A(iV) and U{N) as a function of N, and translating these into parameters 
for a temperature scale, T/v = 7?DA~(^~^)/^, where r/ is an appropriately chosen constant of 
order unity. ^ For the particle-hole symmetric case we take the value of Upp{N) {=Uhh{N)) as 
U{N), and translate this, together with A(A^) and ed{N), into parameters appropriate for a 
temperature scale T^r, and generalize the RPT expression for the impurity susceptibility in 
eq. (12) to finite temperatures. 



Xs{T) = 



Pd(o,r)(i + c/(T)pd(o,r)) 



(23) 



2 
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Fig. 7. The logarithm of a2{h) (stars), the coefficient of the (T/Tk)^ term in the expression (21) for 
the low temperature conductivity (t(T, h), is plotted as a function of the logarithm magnetic field 
h using the renormalized parameters shown in Fig. 6 (symmetric model with U / tt A = 5) 




h/jtA 



Fig. 8. The coefficient of the (T/Tk)^ term (stars) in the expression (21) for the low temperature 
conductivity cr(T, h) is plotted as a function of the magnetic field h over a range Q < h < 1.5Tk, 
using the renormalized parameters shown in Fig. 6 (symmetric model with U/-k^ = 5, Tk/ttA = 
0.00204) 



where Pd{0, T) is the free quasiparticle contribution to the impurity susceptibility given by 

pa{0, T) = ~ p Pd{uj)^^cL: (24) 
where /(w) = l/(e'^/-^ + 1), and Pd{i^) is the free quasiparticle density of states given by 
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Fig. 9. The ratio of Xs(^)/Xs(0) (stars) deduced from eq. (23) using temperature dependent param- 
eters A(r) and U{T), for the symmetric model with U/nA = 5 as a function of the temperature 
T. The full curve is that derived from the Bethe ansatz solution of the s-d model. 

eq. (5). We calculate pd{0,T) in the Kondo regime for U/ttA = 5.0 at values of Tjv, using 
the renormalized parameter A(rAr), with ed{T]\f) = and rj = 1.2, as is used in the NRG 
evaluation of spectral densities on a scale con (see for example^^'^^). We then deduce Xs{T) 
from eq. (23) using U{Tn)- In Fig. 9 we compare the results of this calculation with the 
Bethe ansatz results for the s-d model given in reference. There is excellent agreement with 
the exact Bethe ansatz results over this temperature range. The enhanced value of U (T) at 
very high temperatures, as in the high field case, is qualitatively understandable in terms of 
mean field theory, where Umt{T) = U/{1 — ^7p(i,mf (0, T)), which when substituted into eq. 
(12) gives the mean field value for Xs{T)/{gfiBf = 0.5pd,mf(.^,T)/{l - Upd,mi{0,T)). The 
value of Xs{T)/{gfiB)^ in the extreme high temperature range corresponds to that of the free 
bare model, 1/8T. It can be seen from Figs. 1 and 2 that in the smaller N range there is no 
unique prediction for U{N), so this calculation is not well defined. Nevertheless, the agreement 
with the Bethe ansatz results is remarkable and it does indicate that an RPT expansion with 
temperature dependent parameters might be possible. Such a perturbation expansion, being 
a general technique, would have a potentially wide application to other strongly correlated 
systems. 

6. Conclusions 

We have shown that low energy or strong coupling NRG fixed point of the s-d and Ander- 
son models can be analysed as a renormalized version of the Anderson model. This analysis 
has the advantage that the Fermi liquid aspects, and the 1-1 correspondence of the single 
particle excitations with those of the non-interacting system, is brought out clearly. The three 
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renormalizcd parameters, id, A and U, can be used to specify completely a renormalized 
perturbation expansion, which is applicable on all energy scales. 

With a magnetic field present, renormalized parameters, ed{h), A{h) and U{h) have been 
calculated as a function of the magnetic field strength h. Using these parameters in the 
renormalized perturbation theory, we have derived an expression for the low temperature 
conductivity as a function of magnetic field strength. 
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